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The Somos 5 sequences are a family of sequences defined by a fifth order bilinear 
recurrence relation with constant coefficients. For particular choices of coefficients and 
' initial data, integer sequences arise. By making the connection with a second order 

nonlinear mapping with a first integral, we prove that the two subsequences of odd/even 
index terms each satisfy a Somos 4 (fourth order) recurrence. This leads directly to 
the explicit solution of the initial value problem for the Somos 5 sequences in terms of 
the Weierstrass sigma function for an associated elliptic curve. 
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The subject of this article is the bilinear recurrence relation 

T n +3T n -2 = OLT n+2 T n ^i + (3 T n+ iT n , (1.1) 

where a and f3 are constant parameters, which is known as the Somos 5 recurrence JT2H3JEB]- 
Somos 5 sequences are denned by specifying initial data To, Ti, r 2 , r 3 , r 4 together with the 
values of the coefficients a, j3, and when all of these are suitably chosen then integer sequences 
can result. For example, when 

r = n = r 2 = r 3 = r 4 = 1, a = (3 = I (1.2) 

then the integer sequence 

1, 1, 1, 1, 1, 2, 3, 5, 11, 37, 83, 274, 1217, 6161, 22833, . . . (1.3) 

results. Traditionally the latter is referred to as the Somos 5 sequence, but we will follow 
the terminology of [211122] and refer to this specific sequence as Somos (5). 

In recent work ^3], we have considered fourth order quadratic recurrences of the form 

r„+ 2 T„- 2 = a Tn+lTn-! + (3 (r n ) 2 , (1.4) 

where a and f3 are constant parameters. Sequences that satisfy such fourth order recurrences 
are known as Somos 4 sequences [2111201, and they originally arose in the theory of elliptic 
divisibility sequences [33 EH 123 • in that context, both the parameters a, j3 and the iterates 
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r n are usually integers, in which case the sequences have applications in number theory, as 
they provide a potential source of large prime numbers [HI El- Moreover Somos 4 sequences 
provide a simple example of the Laurent phenomenon: taking the initial data to,ti,T2,T3 
and the parameters a, (3 as variables, all subsequent terms r n for n > 4 in the sequence 
are Laurent polynomials (rather than just rational functions) in these variables. Fomin 
and Zelevinsky have proved that this remarkable "Laurentness" property is shared by a 
variety of other recurrences in one and more dimensions, including the Somos 5 recurrences 
(jl.lj) . and this has applications in combinatorics and commutative algebra (see [11J and 
references). Although Somos sequences are often defined over Q, here we will work over C 
and occasionally indicate when our results restrict to yield rationals or integers. 

Morgan Ward first introduced elliptic divisibility sequences |331 13*4] , a family of antisym- 
metric sequences with a n = — a_ n defined by recurrences of the form 

^n+2^n-2 — ( a 2) 2a n+lfln-l — a l a 3( a n) 2 ; (1-5) 

by considering sequences of rational points nP on an elliptic curve E over Q (see also "I21I2S])- 
To obtain integer sequences of this kind it is required that 

a = 0, a,\ — 1, a 2 ,a 3 ,a 4 GZ with a 2 |a 4 , (1.6) 

and then it turns out that the subsequent terms of the sequence satisfy the remarkable 
divisibility property 

c*n|o m whenever n\m. (1-7) 

These integrality and divisibility property can be said to rest on the fact that the terms of 
an elliptic divisibility sequence satisfy the Hankel determinant relation 



"n+m ""n—m 



(1- 



for all m, n e Z. Starting from the Hankel determinant formula, with initial data of the form 
(jl.fi") . it is easy to prove by induction that all a n for n > are integers with the divisibility 
property (jl.7|) . Although with a$ = there is potentially the problem of dividing by zero, it 
turns out that these sequences are consistently extended to all fiGZ. 

As we explain below, the formula (jl.8j) itself is most easily proved by making use of 
Ward's explicit expression 

«„ = ^. (1.9) 

for the general term of the sequence in terms of the Weierstrass sigma function associated 
with the curve E, and then applying the addition formula 

a(z + k.) a(z — «),,, v . i 

auralr = P{K) - p(z) (L10) 

or its immediate consequence, 

<t(c + d) a{c — d) <r(a + b) a(a — b) 
-a(b + d) a(b - d) a(a + c) a(a - c) (1.11) 
+a(b + c) a(b — c) a(a + d) cr(a — d) = 0, 

which is called the three-term equation (see §20.53 in [H*5])- Ward's elliptic divisibility se- 
quences are just special cases of Somos 4 sequences. The following theorem was proved in 



Theorem 1.1. The general solution of the Somos 4 recurrence relation \l-4\j takes the form 

o-(z + nn) 



r n = AB r 



where k and z are non-zero complex numbers, the constants A and B are given by 



A 



TO 



B 



ct(k)ct(z ) t x 



1.12 



1.13) 



a(z ) ' (t(z + k) t ' 

and o~(z) = cr(;z; g 2 > #3) denotes the Weierstrass sigma function of an associated elliptic curve 



E : y 2 = 4x 3 - g 2 x - g 3 . 



(1.14) 



The values k, zq and the invariants g 2 , g 3 are precisely determined from the initial data 
r o? r i) r 2? T3 and the parameters a, (3, where the latter are given as elliptic functions of k by 



;i.i5) 



The terms of a Somos 4 sequence correspond to the sequence of points Pq + nP on 
the curve E, with zq + nn being the associated sequence of points on the Jacobian. The 
elliptic divisibility sequences arise when P — > 00 (or equivalently z — > with r — > 0). To 
illustrate the effectiveness of such explicit formulae for these sequences in terms of Weierstrass 
functions, we can mention two corollaries of the above theorem. 

Corollary 1.2. The terms of a Somos 4 sequence satisfy the Hankel determinant formula 



7~>n.4-m.7~i 



n+m ' n—m 



^m—l^n 



;i.i6) 



for allm,n G Z, where a m = a (run) / 'a(K) m are the terms of an associated elliptic divisibility 
sequence. 



Proof. Substituting the explicit expressions (jl.9|) and (|1.12|) for a n and r n into (jl.lfcij) . setting 
z = zq + uk and using (|1.10|) yields 



a(K 



,2m 2 +2n 2 



A 2 B 2n a(mK) 2 a(z) 2 



Trt.+m.Ti 



n+m ' n—m 



(a m ) 2 r n+ iT n _i + a m +ia m _i(r n ) 2 \ 



<y{z + vnK)a{z — ran) a(z + n)a(z — k) a((m + l)re)cr((m — 1)k) 



aizYain) 2 



a(K) 2 a(mK) 2 



a(z) 2 a(mK) 2 

= p(mK) - p(z) - p(n) + p(z) + p(n) - p(mK) = 
Corollary 1.3. The terms of a Somos 4 sequence also satisfy 



dm+l'Tn+2 ^m—Y^n 
Q J m+2l~n+l ^mTn—1 



Ul^Tn+m+lTn-m 1 

for all m,n G Z, with a n for n G Z given by JTjJ) as above 



□ 



(1.17) 



Proof. Substituting the formulae (jl.9j) and (jl.l2j) for a n and r n into ()1.17|) . and setting 
z = z + (n + \)k gives 



A 2 B 2n+l 



a lO'2 T n+m+l T n~m ~ &m a m+l T n+2 T n-l + Ojm-lOj m +2% 



= ct(k) o{2k) a[z + (m + 1/2)«J ct^ - (m+ 1/2)kJ 
— a{niK) cr^(m + er^z + 3/t/2^ a^i: — 3k/2^ 
+ cx((m - l)Kj cr((m + 2)/tj er(z + k/2\ a(z - k/2 



by (HHD. 



□ 



Remarks 1.4. These two corollaries of Theorem 1.1 are the main result of [2ES|, where they 
are proved by purely algebraic means i.e. without using explicit analytic expressions such 
as ()1.12j) . Morgan Ward's formula (jl.8|) for elliptic divisibility sequences is clearly just a 
special case of the first corollary. Comparing the statement of Corollary 1.2 when m = 2 
with the expressions for the coefficients of the Somos 4 recurrence in (|1.15J) gives a = (c^) 2 , 
(3 = —aia 3 , so it is clearly required that 

PV) 2 = ^£ P'^f Wk) - P {k)) = (1.18) 

To see why these two identities must hold, note that for any n, a n = <r(n/c) / 'o~(K) n2 is an 
elliptic function of k (see exercise 24 in Chapter 20 of [35j). Thus in the first formula of 
([1.18)1 . both sides are elliptic functions of k whose Laurent expansions around k = have the 
same principal parts; moreover both functions vanish at any one of the half-periods K = Uj, 
j = 1,2,3; hence these two expressions for a are equivalent, and thence the two different 
formulae for (3 are seen to be equivalent by taking z = 2k, in (jl.lUj) . 

The main results of this paper are found in the next section: Theorem 2.7 presents the 
explicit solution of a second order nonlinear mapping associated with Somos 5 sequences, 
while Theorem 2.9 presents the explicit form of the general solution for these sequences. 
In the third section we briefly illustrate the solution of the initial value problem for such 
sequences, by applying the theorem to the specific example of the Somos (5) sequence (jl.3|) . 
which is sequence A006721 in Sloane's catalogue 



2 General solution of Somos 5 

The Somos 4 recurrence has a simple gauge invariance property: if r n satisfies (11.4)1 then for 
any non-zero constants A, B it is clear that the transformation 

r n ^f n = AB n Tn} neZ (2.1) 

yields another solution of the same fourth order recurrence. Thus it is natural to introduce 
the gauge-invariant combination 

r _ T n+l T n-l /„ 9 n 

in ~ {2.2) 

which satisfies the second order nonlinear mapping 



(written in the more symmetrical form (|2.4j) below). The following result in ^3] was the 
basis for the proof of Theorem 1.1. 

Proposition 2.1. The second order nonlinear mapping 

fn-l(fn) 2 fn + l = af n + P (2.4) 

has the conserved quantity 

J = fn-lfn + a + + ~r^~~r- (2-5) 



fn J fn—lfn 

Its general solution has the form 

fn = p(k) - p(z + rut), (2.6) 

where p(z) = p(z;g2,gz) is the Weierstrass elliptic function associated with an elliptic curve 

\l-14\j with invariants g<i, #3. The parameters z , k e C and the two invariants are deter- 
mined from a, (3 and the initial data fo, f± according to the formulae 

X^^-(L-/3\ g 2 = 12X 2 -2J, g 3 = 4X 3 -g 2 X-a, (2.7) 



3a V 4 

together with the elliptic integrals 



: x - fo dx r x dx 

z = ± — « = ± / — . (2.1 



y Joo y 

Remarks 2.2. The above is a slight reformulation of Proposition 2.2 in ^Hj, where the quan- 
tity A = p(k) was taken as the basic conserved quantity rather than the simpler quantity 
J = p"(k)\ the latter is the same as the "translation invariant" found in Swart's thesis |2*§] . 
There is an ambiguity in the sign of the integrals ()2.8|) due to the elliptic involution y — > —y, 
but the solution ()2.6|) is clearly invariant under Zq — > — Zq, K — > — K due to the fact that 
p(z) is an even function; the relative sign of zq and k is fixed by the consistency require- 
ment p\zq)p'(k) = (/o) 2 (/i — The mapping (|2.4jl . or equivalently (|2.Hjl . is of the type 
studied by Quispel, Roberts and Thompson [2H] (see also Pfl HU [T31 EO] and references). The 
solution ()2.6|) of (|2.4|) corresponds to the solution ()1.12|) of the Somos 4 recurrence, modulo 
the two-parameter abelian group of gauge transformations (|2.1j) . 

Proposition 2.3. By means of the substitution the Somos 5 recurrence M.l)) corre- 

sponds to the third order nonlinear mapping 

fn-l(fn) 2 (fn + l) 2 fn+2 = &f n f n+1 + 0, (2.9) 

which has two independent conserved quantities I , J given by 

I = fn-lfnfn+X + « ( ~F~ + T + tM + f J f . ( 2 -10) 
\Jn— 1 In Jn+l / In—llnln+l 

J = fn-lfn + /n/n+1 + 6t ( — + — ) + ~? ( f \2 f ' i 2 - 11 ) 

\/n— 1/n Inln+lJ In—\\In) Jn+l 

Any solution of the Somos 4 recurrence \l-4\ l also satisfies a Somos 5 recurrence M.l)) with 
parameters 

& = -/3, f3 = a 2 + (3J, (2.12) 

where J is as in Met. 5)) . and values of the conserved quantities for the associated solution of 
i TOj) are 

J=J, I = 2a. (2.13) 



Proof. Letting I n denote the expression on the right hand side of (|2.10j) . and letting I n+ % 
be the same expression with n — > n + 1, we have immediately that I n+ \ — I n — whenever 
()2.9j) holds, and similarly J in (|2.11j) is conserved by the third order nonlinear mapping. 
Now for any sequence r n , n G Z that is a solution of a Somos 4 recurrence (jl.4j) . there is a 
corresponding solution of ()2.4|) given by ()2.2|) . By considering the expression 

£n '■= fn-l(fn) 2 (fn+l) 2 fn+2 ~ &fnfn+l ~ $, 

and then using (j2.4j) , first with n — > n + 1 to eliminate / n+2 , and then once more to eliminate 
f n+1 , we find 



£ n = a 2 -/3 + /3/ n /„_!+/3 



/n fnfn—lj fn—1 



a 2 + (3 J - (3 - (/3 + a) ( -7^- + ' 



K fn—1 fn~lfn, 

(where we have used (J2.5J) at the last step) and hence £ n = for all n only when the two 
conditions (j2.12|) hold. In that case, f n satisfies (|2.9J) with these values of a, /?, which is 
equivalent to r n being a solution of The values (|2.13|) follow from similar calculations. 

□ 

Remark 2.4. Clearly the special case m = 2 of Corollary 1.2 above means that every Somos 
4 is a Somos 5 with a = a 3 , (3 = — a 4 /a 2 , so using ()1.9j) to compare with (|1.15j) and noting 
that J = p"{k) it is clear that the equation a = a 3 = —(3 is the same as the first equation 
of (jl,18|) . while the formula for (3 in (J2.12|) yields the interesting identity 

P \kY + p"(k) P \k) 2 ( P (2k) - P (k)) = - ^w\i2 - ( 2 - 14 ) 

The latter should be very tedious to prove by comparing poles and zeros, since on each side 
the leading term in the Laurent expansion at n = is — 2/k 12 . 

The formula (jl,12j) is obviously not general enough to capture all of the Somos 5 se- 
quences: not every Somos 5 is a Somos 4! The expression (jl,12j) only depends on the six 
parameters A, B, z , k, g 2) whereas to specify a generic Somos 5 sequence requires seven 
parameters, corresponding to the five initial data r , 7~i, r 2 , r 3 , r 4 plus the two coefficients 
a, (3. One way to interpret the extra degree of freedom in Somos 5 sequences is by noting 
that they have a three-parameter group of gauge transformations given by 

T 2n — > r 2n = A + B 2n r 2n , T 2n+1 — ► f 2n+ i = i_ B 2n+1 r 2n+1 , n G Z. (2.15) 

In other words, compared with Somos 4 there is an extra freedom to rescale the even/odd 
terms differently while leaving the recurrence (jl.ip unchanged. The quantity defined by 
()2.2|) is not invariant under this larger group of transformations, and thus it is natural to 
consider something that is, namely 

K = fn+ifn = Tn+2Tn ~ 1 , (2.16) 

Before we consider the recurrence satisfied by h n , it is convenient to state the following 
general result on mappings of the plane associated with biquadratic curves, which have been 
studied in detail recently by Iatrou and Roberts [T^ ITo] , Bastien and Rogalski lj and also 
Tsuda |30| . 



Proposition 2.5. The quartic curve E defined by 



B(X, Y) = eX 2 Y 2 + dXY(X + Y) + c{X 2 + Y 2 ) + b(X + Y) + a- KXY = (2.17) 
has genus one. The mapping of the plane 

(u n -i,u n ) i-> (u n ,u n+1 ) (2.18) 

defined by 

a + bu n + cu 2 n 

Un+lUn-l = — — } : 2 l 2 ' 19 ) 

c + du n + eu n 

corresponds to the composition of two involutions on E, and has the conserved quantity 
K = eu n -iu n + d(u n -i + u n ) + c + — ^ +b[ + — + • (2.20) 

V U n U n -1 J \U n -l U n J U n -\U n 

The iterates of correspond to the the sequence of points 

P n = (w , Mi) + nP E E (2.21) 
for some P G E. Furthermore, upon uniformizing E, the coordinates u n can be written as 

u n = f(z + UK) (2.22) 

in terms of an elliptic function f and certain parameters Zq, k G C. 

Proof. To see that the curve E is elliptic, it is sufficient to observe that by the transformation 
1Z = XY, S = X + Y, this symmetric biquadratic is a ramified double cover of the rational 
curve 

eTZ 2 + dTZS + c(S 2 - 211) + bS + a - KK = 

with branching at four points (X, X) given by the roots of the quartic B(X, X) = 0, and 
so by the Hurwitz formula the genus is one. The symmetric biquadratic (J2.17j) admits the 
involutions 

iq: (X,Y)^(Y,X), F 2 : (X, Y) i-> (X, Y^), 

where y — Y^ is the unique solution to the quadratic equation B(X, y) = apart from y — Y. 
Clearly the composition G — F 2 o F\ of these two involutions yields the mapping of the plane 
(|2.18|) defined by (|2.19|) . and it is easy to see that the quantity (|2.2Uj) is conserved by this 
mapping, so that K = K = constant and the points (u n -i, u n ) lie on the level set (|2.17j) for 
all n (see e.g. PP for more details of this construction). G is thus a birational map E — > E 
composed of two involutions Fj, j = 1, 2. We claim that for all P G E each Fj : 2Cj — P 
for some point Cj G E, for j — 1, 2 respectively, and hence G is a translation by the point 
P = 2C 2 - 2Ci, so G : P i-> P + P as is required for TTT^ to hold. 

To see why the claim is true, note that each involution Fj, j — 1, 2 is a birational map 
E — ► E with four fixed points. Consider first F±, noting that any morphism of elliptic curves 
is a combination of a translation by a point and an isogeny (see p.75). Since F± is 
birational we require that G Aut(E), and because it is an involution the only possibilities 
are = [±1]. If = [1] then we must have F\ : P k P + Q where f2 is a non-identity 
element of the 2-torsion subgroup (of order four), but then there are no fixed points. Thus 
it must be the case that = [—1] and so F\ : P i— > 2C*i — P where the four possible points 
C\ + Q (with being 2-torsion) account for the fixed points, which for F\ correspond to the 
roots of B(X,X) = 0. The same argument applies to F 2 . 



To obtain the expression (|2.22j) . observe that because E has genus one it can be uni- 
formized using elliptic functions; in [3S] this result is attributed to Clebsch. Thus we can 
write (X,Y) = (g(z), f(z)) where g,f are rational functions of the Weierstrass functions 
p(z), p'{z) associated with a curve E in Weierstrass canonical form (jl.l4|) with invariants 
g 2 , <?3, and z G J&c(E). Now the map defined by (j2.19|) corresponds to addition of a point, 
which is equivalent to a shift z \— > z+^on the Jacobian. Then we see that for some z n we have 
{u n -i,u n ) = {g{z n ) , f(z n ) ) and (u n ,u n+1 ) = (g(z n +i) , f(z n +i) ) = {g{z n + «) , f(z n + «) ), 
whence it follows that g(z) = f(z — n) and w n = g(zo + (w+ = /(-2o + ^ K ) as required. □ 

Remark 2.6. In essence, the above result for symmetric biquadratic curves was already known 
to Euler (see Theorem 9 in [32J). A thorough treatment of general integrable maps associated 
with elliptic curves over an arbitrary field can be found in the recent work For our 

purposes here, we need to construct the explicit uniformization for particular families of 
curves of the above type, and the exact form of the corresponding expressions in terms of 
Weierstsrass functions changes when one or another of the parameters in (j2.17|) degenerates 
to zero. For example, we see that for the mapping (|2.4|) associated with Somos 4 sequences, 
the invariant curve defined by (|2.5jl in Proposition 2.1 corresponds to the degenerate case of 
(j2.17|) obtained by setting e = 1, d = c = 0. In the proof of the analogous result for Somos 
5 sequences, we shall see that the relevant mapping corresponds to the case when (J2.17I) 
degenerates to a cubic curve by setting e = c = 0, cf. (|2.24j) and (|2.34|) below. 

Theorem 2.7. Given a Somos 5 sequence satisfying M.l\) . the corresponding quantity h n 
defined by \2. 1 6]) satisfies the second order nonlinear mapping 

h n -ih n h n+1 = ah n + ft, (2.23) 

which has the conserved quantity 

j = hn _ x + hn + & f _L + _L ) + -JL-. (2.24) 
\n n -i n n J h n -\h n 

The general solution of this map can be written as 

h _ or(*o + (w + 2)/c) <r(z + (w - _ p'(k) f p'(zo + nK)- p'(k) \ | p"(n) 
o-(k) 4 o-(zo + nn) cr(z + (n + 1)k) 2 \ p(z Q + nn) - p(n) J 2 

where a(z) = a{z]g2,g^) and p(z) = p{z]g 2 ,g^) are Weierstrass functions associated with 
an elliptic curve 

E : y 2 = 4x 3 - g 2 x - g 3 . (2.26) 

The invariants of the curve are determined from the parameters a, (3 together with the initial 
data ho, hi for the map \2. 2$) according to the formulae 



fr=(P + aJ) 1 /\ \= — l±- + ck\, (2.27) 




g 2 = 12A 2 - 2 J, g 3 = 4A 3 - g 2 \ - fx 2 . (2.28) 
The parameters z , kgC are then determined by the elliptic integrals 

f X0 dx f~ x dx , n , 

z = ± — , k = ± / — . (2.29) 



00 



y Joe y 



on the curve E, where 

x = X + - ~, h -i = —TT — (2-30) 

h-i + h - J h hi 

and the relative sign of k and z$ is fixed by the requiring the consistency of 

£ = p'(K,)p'(z ) = (x - A)(/i_i - ho). (2.31) 

Proof. Letting J n denote the right hand side of (|2.24|) . it is clear that 

J n +i - J n = ~t~~~t f"" 1 ' \ h n +\h n h n _ x - oth n - (3), (2.32) 

"n— l"n" , n— 1 

so assuming that h n+ \ ^ h n -i for all n then J is conserved if and only if h n satisfies ()2.23j) . 
Alternatively note that rewriting (|2.11|) in terms of h n gives precisely (j2.24j) . and rearranging 
the latter expression for J yields 

(/i n _! + K-J) {h n -xh n + a)+(3 + aJ = 0, (2.33) 

which implies that for any sequence h n , n G Z satisfying ()2.23|) . (X, Y) = (h n _i,h n ) is a 
point on the curve 

(X + Y- J)(XY + a) +{3 + aJ = (2.34) 

for all n. 

We claim that this elliptic curve is birationally equivalent to (j2.26J) via the transformation 
x = E(y±i)+l Y = 4(^£]+l (2.35) 



2\x-\J 2' 2\ X -\J 2' 

where the parameters A, jx and invariants g^, gz are found from the coefficients of the curve 
(J2.34)) using the formulae (J2.27j) and (|2.28|) . To verify this, it is instructive to use the analytic 
formulae (J2.25j) which provide a uniformization of (|2.34|) . First note that the two formulae 
for h n in ()2.25|) are easily seen to be equivalent by noting that both expressions are elliptic 
functions of zq with simple zeros at zq = —(n + 2)n, — (n— 1)k and simple poles at zq = —nn, 
— (n + 1)k, with the same residues at these poles. Similarly, shifting n — >• n — 1 in the first 
(sigma function) expression in 1)2.25)1 gives 

_ Q-{zq + {n + 1)k) a{z + {n - 2)k) _ p'(n) f p'(z + nK) + p'(k) \ p"(k) 
a(K) 4 a(z + (n - 1)k) a(z + nn) 2 \ p(z + uk) - p{K) J 2 

where once again these two formulae for h n -\ are seen to be equivalent by looking at their 
poles and zeros. Then we can make the identifications 

A = P (k), jx = p'(K) = + &jy/\ J=p"(K), & = -p\k) 2 (p{2k)-p{k)), (2.37) 

which are consistent with ()2.27j) and ()2.28|) . i.e. (A, fi) = (p(n), p'{k)) is a point on the curve 
(12~2T)1) . 

Taking the sum of the right hand (p function) expressions in ()2.25|) and ()2.36)) gives 

K-! + K = p"{k) + —, (2.38) 

p{Z + UK) - p{K) 



which by ()2.37|) is equivalent to X + Y = J + fi 2 /(x — A) with the identification (x,y) = 
(p(zq + uk), p'{zq + nn)). On the other hand, taking the product of the left hand (sigma 
function) expressions in (|2.25jl and (j2.3(ij) gives 

a(z + (n + 2)k) a(z + (n - 2)k) , 2 / - . - a , . 

/i n _i/i„ = — = p («) I p(2k) - p(z + rwc) J , (2.39) 

where we have used (11.10)1 and the first identity in ()1.18j) . With the same identifications 
as before, this product identity implies that we have XY + at = — fi 2 (x — A), and hence 
(X + Y — J)(XY + a) = = —(/3 + a J) which means that (X, Y) = (/i n -i, h n ) lies on 
the curve (J2.34J1 for any n, as claimed. In turn this means that the sequence {h n \n G Z} 
defined by ()2.25j) has J = p"(k) as a conserved quantity, and hence by ()2.32|) it must satisfy 
the second order map ([2.23)1 . 

To see that (|2.25j) represents the unique solution of the initial value problem for the map 
()2.23j) . observe that it depends on the four parameters g 2 ,g3,zo,K. Up to the involution 
Zq — > — Zq, K — > — K (which does not change the solution), all of the latter are uniquely 
determined from the initial data ho, hi and parameters a, j3 by first constructing the first 
integral J = Jq (i.e. setting n = in (|2.24jl ) and then using the formulae (|2.27)l and ()2.28|1 
to find the invariants. Finally the base point zq and the shift k on the Jacobian Jac(E') 
of the curve ()2.26|) are obtained by inverting xo = p(zo), A = p(k) to yield the elliptic 
integrals (|2.29|) . where Xq in (|2.30J) is found by rearranging (|2.38|) for n = 0. The consistency 
requirements ()2.31|) follow from fixing the branch of fourth root taken for \x in (|2.27j) and 
ensuring that the value of p'(zq) agrees with this in the formulae for Hq and h^\ (that is, the 
p function expressions in (|2.25j) and (|2.36p when n = 0). □ 

Proposition 2.8. For any solution r n , n G Z of a Somos 5 recurrence M.l\) . the subsequences 
T n = r 2n and t* = T2n+\ } °f respectively even and odd index terms, both satisfy the Somos 4 
recurrence 

< +2 <-2 = + P K)\ (2-40) 

where 

a* = (3 2 , p* = a(2(3 2 + a(3J + a 3 ), (2.41) 
with J being the associated integral \2. 2$ of the map 
Proof. Looking at the formula (|2.39|) we see that 

h n h n _ x = / n+ i(/ n ) 2 /n-i = r " +2T "~ 2 = p'(k) 2 (p(2k) - p(z + uk)), 
which is almost identical to the formula ()2.6|) for the solution of the map (|2.4|) . except for 



the scale factor /i 2 = p'(n) 2 = y (3 + aJ. Using the scaling property 

p(z; 92, 93) = jl~ 2 p(z/fL; fi A g 2 , /i 6 ^), o~(z; g 2 , 93) = ficr(z/fl; fi 4 g 2 , /i 6 ^ 3 ), 

of the Weierstrass functions, we can rewrite everything in terms of rescaled sigma and p 
functions with invariants g\ = fi 4 g 2 , g% = fi 6 g3, setting v — k/Jjl and u = Zq/Jjl to obtain the 
canonical form 

( n 2 = P( 2v ) - p( u o + nv) = — • (2.42) 

It follows immediately from Proposition 2.1 that both F* = r 2n+2 r 2n _ 2 / (r 2n ) 2 and F* = 
T 2n+3 T 2n-i/ ( T 2n+i) 2 ar e solutions of the nonlinear map 



F n+l {F n ) 2 F n _ 1 = a*F n + (3\ 



with 

a* = p'(2v)\ (3* = p'{2vf (p{Av) - p{2v)) . (2.43) 

The required result follows immediately, upon noting that the duplication formula for the 
Weierstrass p function can be used in ()2.43ft to rewrite p(4t>) in terms of p(2v ) and its 
derivatives, and similarly for p(2v) in terms of p(v). After further tedious computations, 
scaling back to the original Weierstrass functions a = p(n), fi = p'(k), J = p"(n) with 
invariants g 2 , #3, we arrive at the expressions (j2.41|) . □ 

Theorem 2.9. The general solution of the Somos 5 recurrence takes the alternating form 

^ = A + K^±^, r w = A. B» ^ + < 9 2 lt 1)V \ (2.44) 

where here c(u) = a(u; g 2 , g%) denotes the Weierstrass sigma function associated with the 
elliptic curve 

E* : y 2 = 4x 3 - g\x - g* 3 , (2.45) 

where Uq, v and the invariants are given in terms of the parameters \2.21 ), \2.2fy) in Theorem 
2.7 by 

92=P l 92, 93 =^93, U = Z /fl, V = K/ji, 

and the pref actors are 

A + = -r-v A - = (\ \ ^ 

a(u ) a(u + v) 

B ^ = a{2v)a{u )T 2 fl = q-(2v)o-(u + v)r 3 ^ 
cr(u + 2v)t ' a(uo + 3v)ti 

The solution \2.44l) depends on seven independent parameters g 2 , g%, u , v, A + , A_ and B + , 
with B± defined in \2.47\) being related by 

|+ = -a(2v). (2.48) 

Proof. Starting from the result of Proposition 2.8, the formula (|2.42j) for the ratios of the r n 
can be used to solve for the even and odd index subsequences separately, since these are both 
Somos 4 sequences. In that case, Theorem 1.1 applies separately to the two subsequences 
r 2n and r 2n+ i, which yields (12.44ft - the proof for each subsequence is the same as in [HI] . 
However, the solution can depend on at most 7 parameters, as is confirmed by the relation 
(|2.48ft between B±. To prove the latter relationship, it suffices to rewrite the ratios of sigmas 
and r n in terms of h\, and see that most terms cancel, taking care to scale appropriately 
with 

fi> = p'(^92,93), /i 4 = p'{v;gl,gl). (2.49) 

□ 

Remark 2.10. The rescaled curve (|2.45ft is more convenient when working in fields other than 
C, since then one does not need to define the curve over the field extension by the fourth 
root p. = (f3 + aJ) 1/A . 

Corollary 2.11. The solution of the third order nonlinear map pO ) can be written in the 
alternating form 

f 2 = /„ ( P(^) ~ p(u + 2nv) \ +i = h( PW"^^^ 11 )^ ( 2 50) 

2n ° V p( v ) - p( u o) J ' 2n+1 1 V p(v) ~ p( u o + v ) / 

where p{u) = p(u; g 2 , g£) is the Weierstrass p function associated with the curve \2.4ty , with 
the parameters uq, v as given in Theorem 2.9. 



Proof. This follows by substituting (|2.44jl into (|2.2j) and rewriting suitable combinations of 
the prefactors A±, B± in terms of fo and f\. □ 

Corollary 2.12. The terms r n , n e Z of a Somos 5 sequence satisfy the Hankel determinant 
formula 



a l a 2 T n+m+l T n- 



:2.5i) 

for all m,n G Z ; where the quantities a m , which form an associated elliptic divisibility 
sequence, are given by 

aOmn) „ 
flm = , w , m e Z 

m terms of the sigma function o~(z) = o~(z; g2, #3) /or £/ie curve \2.2b)) in Theorem 2.7. 
Proof. Using the scaling by ft as in (|2.49|) . the formulae (|2.44|) can be rewritten as 

r 2n = A + B + ——2 , T 2n+ i = A- B_ ——2 , (2.52) 

a(2K) n o\Ik) u 

in terms of the sigma function for the original curve (|2.2fjjl as in Theorem 2.7. In terms 
of this sigma function it is easy to show, upon rescaling ()2.48|) by ft — p'(k; g 2 , gs) = 
— a (2k; g 2 , g3)/o-(n; g 2 , O3) 4 (compare with the first formula of (|1.18|l ). that 

B+/B- = o-(K;g 2 ,g 3 ) 4 . 

For verification of (|2.5ip . the symmetry m — > — m — 1 means that it is sufficient to check 
only the cases when m, n are either both odd or both even, to take the different expressions 
()2.52j) for r n with odd/even n into account. Thereafter the proof is the same as for Corollary 
1.3. □ 

Corollary 2.13. The terms of a Somos 5 sequence, satisfying the recurrence M.l)) . have the 
leading order asymptotic behaviour 

log|r n | ~ file - -\og\a(2v)\] n 2 , n -»• 00 (2.53) 



2cji J 4 

for v G UiM., where it)\ is a half-period and r)\ = C(^i); with the invariants g 2 , o| being as in 
Theorem 2.9. 

Proof. The sigma function is given by the product formula 

. , 2u) x (r} X z 2 \ . ( ttz \ tt f 1 2q 2n cos(nz/uj 1 ) + q An ] 

° {z) = - exp (2^) sm (^J n { (T^y } 



with q = exp(Triu 2 /u}i) (see §20.421 and also §21.43 in for the corresponding theta 
function expression). The result then follows upon taking the logarithm of (j2.44j) . The 
corrections to (|2.53jl at 0(n) can be calculated similarly, although these have a different 
form for even/odd n. □ 



3 Example: Somos (5) 



To illustrate the above results, we consider the example of the Somos (5) sequence (jl.3j) . 
The initial data and parameters as in (jl.2|) give ho — 2, h\ — 1 for the second order map 
()2.23|) . whence h-i = 3/2 with a — 1 = j3, so applying the formulae in Theorem 2.5 we have 

,7=5, jfj, = b ' , A = p, x 



12V6' U 12V6' 

whence g 2 = 121/72, g 3 = — 845/(1296\/6). Rescaling by appropriate factors of yEt, as in 
Theorem 2.9, we arrive at the curve 

2,3 121 845 , 121 „ 845 

75 : u = 4x xH , o 9 = , qZ — , 3.1 

y 12 216 y2 12 ' ya 216 V ; 

which has the j-invariant 

1728(^) 3 1771561 
3 ~ (g*f - 27(g*) 2 ~ 612 

Hence the Somos (5) sequence (jl.3|) corresponds to the sequence of points 

(x*, y*) + n(\*, fi*) = (-19/12, 2) + n(29/12, 6) = (5/12, 0) + (n - 2) (29/12, 6) (3.2) 

on the rescaled curve ()3.1|) . 

For comparison with results of other authors, we note that van der Poorten JHj has 
obtained the Somos (5) sequence from the minimal model V 2 + UV + bV = U 3 + 7U 2 + 
12U, while Zagier [HE] found the related curve y 2 + xy = x(x — l)(x + 2); both of these 
curves have the same j-invariant as (|3. ljl . and hence they are birationally equivalent to each 
other. Zagier's discussion makes use of the inspired substitution ()2.16|) for the Somos (5) 
sequence (jl.3J) . although the presentation of the associated curve in the form ()2.34|) contains a 
typographical error in jHEj. Other approaches to the sequence (jl.3|) can be found on postings 
to Propp's "robbins" forum [22] - see in particular the references to unpublished work of 
Elkies in (H| , where the sequence (11.3)1 is related to the problem of Heron triangles with two 
rational medians. 

With the curve ()3.1|) and the sequence of points ()3.2|) we can use Theorem 2.9 to calculate 
the general term of the Somos 5 recurrence in the form ()2.44|) with 

u = 0.163392411 + u 2 , v = -0.672679183, 

where the real and imaginary half-periods are 

ut = 1.181965956, uj 2 = 0.973928783 i 

respectively (doing all the elliptic integrals to 9 decimal places with MAPLE). The result 
of Corollary 2.13 then implies that the terms of the sequence (|1.3J) have the leading order 
asymptotics 

log r n ~ 0.071626946 n 2 , n -> oo. 



4 Conclusions 

The theorems presented in ^3] and above provide a simple algorithmic approach to de- 
termining elliptic curves associated with Somos 4 and Somos 5 sequences, and finding the 
complete solution of the associated initial value problems. Our approach is based on making 



the connection with explicitly solvable maps with conserved quantities i.e. (|2.4|) and (|2.23jl . 
A symplectic map with sufficiently many conserved quantities is called integrable [SI . As 
they stand the maps presented here are not symplectic, but by making a change of variables, 
a suitable symplectic structure was given in [T3|. 

The results presented above naturally lead to the question of how to solve Somos (2iV + 2) 
recurrences with N + 2 terms, of the form 

N 

T n+N+l T n-N-l = ^ ] &j T n+j r n-j, 
3=0 

with constant coefficients aj, and similarly for Somos (2N + 3). By using the continued 
fraction expansion of the square root of a sextic, van der Poorten [21] has generated a certain 
class of Somos 6 sequences associated with genus 2 curves. With a different approach [2j 
based on the addition formulae for Kleinian sigma functions (see [Hj and references) , we have 
constructed Somos 8 sequences derived from genus 2 curves, taking a different (quintic) affine 
model compared with van der Poorten. The recurrences in J2] generalise the genus 2 case 
of Cantor's hyperelliptic division polynomials 7] (for analytic formulae, see (E|), and they 
are related to a family of integrable symplectic maps (discrete Henon-Heiles systems). We 
should also remark that Buchstaber and Krichever have derived bilinear addition formulae 
for Riemann theta functions 0], which have exactly N + 2 terms in genus N. Whether this 
result can lead to an effective solution to the initial value problem for these higher order 
recurrences is the subject of further investigation. It would also be interesting to see how 
some of our results on Somos 4 and 5 sequences could be modified in the p-adic setting j27J . 
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